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• Abstract 

The mass of the lightest CP-even Higgs boson of the minimal supersymmetric extension 
^ I of the Standard Model (MSSM) has previously been computed including 0{aas) two- 

d \ loop contributions by an on-shell diagrammatic method, while approximate analytic 

results have also been obtained via renormalization-group-improved effective potential 
and effective field theory techniques. Initial comparisons of the corresponding two-loop 
results revealed an apparent discrepancy between terms that depend logarithmically 
on the supersymmetry-breaking scale, and different dependences of the non-logarithmic 
terms on the squark mixing parameter, Xf. In this paper, we determine the origin of 
these differences as a consequence of different renormalization schemes in which both 
calculations are performed. By re-expressing the on-shell result in terms of MS pa- 
rameters, the logarithmic two-loop contributions obtained by the different approaches 
are shown to coincide. The remaining difference, arising from genuine non-logarithmic 
two-loop contributions, is identified, and its effect on the maximal value of the lightest 
CP-even Higgs boson mass is discussed. Finally, we show that in a simple analytic 
approximation to the Higgs mass, the leading two-loop radiative corrections can be ab- 
sorbed to a large extent into an effective one-loop expression by evaluating the running 
top quark mass at appropriately chosen energy scales. 



1 Introduction 



In the minimal supersymmetric extension of the Standard Model (MSSM), the mass 
of the lightest CP-even Higgs boson, m^, is calculable as a function of the MSSM 
parameters. At tree-level, rrih is a function of the CP-odd Higgs boson mass, m^, and 
the ratio of vacuum expectation values, tan (3. Moreover, the tree-level value of rrih is 
bounded by rrih < ^z \ cos 2/5| , which is on the verge of being ruled out by the LEP Higgs 
search |Q. When radiative corrections are taken into account, ruh depends in addition on 
the MSSM parameters that enter via virtual loops. The radiatively corrected value of 
depends most sensitively on the parameters of the top-squark (stop) sector: the average 
squared-mass of the two stops, M|, and the off-diagonal stop squared-mass parameter, 
rritXt. The stop mixing parameter is Xt ^ At — /icot/3, where At is the coefficient 
of the soft-supersymmetry-breaking stop-Higgs boson tri-linear interaction term and /x 
is the supersymmetric Higgs mass parameter. The radiatively-corrected value of m\ 
is enhanced by a factor of Gpfnj and grows logarithmically as Ms increases [0]. In 
particular, the upper bound for ruh (which is achieved when ^ mz and tan j3 ^ 1) 
is significantly increased beyond its tree-level upper bound of mz- 

The complete one-loop diagrammatic computation of rrih has been carried out in 
refs. [^^. However, for Ms ^ rrit, the logarithmically enhanced terms are significant 
(in particular, the most significant logarithmic terms are those that are enhanced by 
the Gprnt pre-factor noted above), in which case leading- logarithmic corrections from 
higher-loop contributions must be included. These terms can be summed via renormal- 
ization group techniques. The result of these corrections is to reduce the one-loop upper 
bound on rrih. For Ms < 0{1 TeV), it is found that rrih ^ 125 GeV, where the maximum 
is reached at large and tan P when Ms is maximal and Xt ~ ±\/6Ms (the so-called 
"maximal- mixing" value for stop mixing). However, sub- leading two-loop corrections 
may not be negligible, and a more complete two-loop computation is required. 

The full diagrammatic calculations lead to very complicated expressions for the ra- 
diatively corrected value of rrih. Effective potential and effective field theory techniques 
have been developed which can extract the dominant contributions to the Higgs mass 
radiative corrections (when Ms is large), resulting in a simpler analytic expression for 
rrih. These methods also provide a natural setting for renormalization group improve- 
ment. Although the exact solution of the renormalization group equations (RGEs) must 
be obtained numerically, the iterative solution of the RGEs can easily yield simple ana- 
lytic expressions for the one-loop and two-loop leading logarithmic contributions to rrih. 
These leading logarithms can also be obtained by expanding the complete diagrammatic 
results in the limit of Ms ^ rrit, and this serves as an important check of the various 
computations. 

The effective potential method ^ provides an important tool for evaluating the 
Higgs mass beyond the tree-level. It can be used to provide a short-cut for the cal- 
culation of certain combinations of Higgs boson two-point functions that arise in the 
diagrammatic computation. The effective field theory (EFT) approach ||^ provides a 
powerful method for isolating the leading terms of the Higgs mass radiative corrections 
when Ms rrit. In this formalism, one matches the full supersymmetric theory above 
Ms to an effective Standard Model with supersymmetric particles decoupled below M5. 
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The Standard Model couplings in the MS scheme are fixed at Ms by supersymmetric 
matching conditions. Standard Model RGEs are then used to evolve these couplings 
down to the electroweak scale (either rrit or mz)- While the stops are decoupled at 
Ms, the stop mixing (so-called "threshold") effects are incorporated by modifying the 
matching conditions at Ms- With the use of effective potential techniques, the EFT 
formalism and the iteration of the RGEs to two-loops, the leading contributions to the 
radiatively-corrected Higgs mass was obtained in analytic form in refs. [p|-pJ]l These 
results included the full one-loop leading logarithmic corrections and one-loop leading 
squark-mixing threshold corrections, the two-loop leading double-logarithmic correc- 
tions and the two-loop leading logarithmic squark-mixing threshold corrections up to 
0{hlas) and 0{hf), where ht is the Higgs-top quark Yukawa coupling. 

In order to extend the above results, genuine two-loop computations are required. 
The first two- loop diagrammatic computation was performed in ref. [|T2| in the limit of 
mA ^ rnz and tan (3 ^ 1 (where rrih attains its maximal bound), where only terms 
of 0{h^) and 0{hfas) were evaluated, and all squark mixing effects were neglected. 
More recently, a more complete two-loop diagrammatic computation of the dominant 
contributions at 0{aas) to the neutral CP-even Higgs boson masses has been per- 



formed |13-15]. This result was obtained for arbitrary values of tua, tan/3 and the 



stop-mixing parameter Xf. This two-loop result, which had been obtained first in the 



on-shell scheme, was subsequently combined [|14|,p!5| with the complete diagrammatic 
one-loop on-shell result of ref. and the leading two-loop Yukawa corrections of 0{hf) 
obtained by the EFT approach P|-pT|. The resulting two-loop expression was then 
expressed in terms of the top-quark mass in the MS scheme. By comparing the final 
expression with the results obtained in refs. it was shown that the upper bound 

on the lightest Higgs mass was shifted upwards by up to 5 GeV, an effect that is more 
pronounced in the low tan j3 region. 

Besides the shift in the upper bound of rrih, apparent deviations between the ex- 
plicit diagrammatic two-loop calculation and the results of the EFT computation were 
observed in the dependence of rrih on the stop- mixing parameter Xt. While the value of 
Xt that maximizes the lightest CP-even Higgs mass is {Xt)ma.x — ±VQ Ms ~ ±2.4 



in the results of refs. p|-p]T|, the corresponding on-shell two-loop diagrammatic compu- 
tation found a maximal value for rrih at (Xf)max ~ SM^.Q Moreover, in the results of 
refs. rrih is symmetric under Xt —Xt and has a (local) minimum at Xt = 0. 

In contrast, the two- loop diagrammatic computation yields rrih values for positive and 
negative Xt that differ significantly from each other and the local minimum in rrih is 



shifted slightly away from Xt = |T^. A similar conclusion was reached in ref. 



which used an effective potential calculation to extend the results of ref. [|12[ to the case 
of non-zero stop mixing. 

A closer comparison of results of the EFT computation of the radiatively-corrected 
Higgs mass and the two-loop diagrammatic computation at first sight revealed a sur- 
prising discrepancy. Namely, the two-loop leading logarithmic squark-mixing threshold 
corrections at Olh'fas) of the former do not appear to match the results of the latter. 
In this paper, we shall show that this apparent discrepancy in the leading- logarithmic 



-•^A local maximum for is also found for Xt ~ —2AIg, although the corresponding value of rrih 
at Xt ~ +2Ms is significantly larger 
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contributions is caused by the different renormalization schemes employed in the two 
approaches. While the original two-loop diagrammatic computations of ref. |]13[ were 
performed in an on-shell scheme, the results of the EFT approach are most naturally 
carried out in the MS scheme. In comparing results obtained within different renor- 
malization schemes in terms of the (not directly observable) parameters Xt and Ms, 
one has to take into account the fact that these parameters are renormalization-scheme 
dependent. The effect of this scheme dependence first enters into the calculation of 
at the two-loop level. In order to allow a detailed comparison between the results of the 
different approaches we derive relations between these parameters in the two different 
schemes. We apply these relations to re-express the diagrammatic on-shell result in 
terms of MS parameters. In this way we show that the leading logarithmic two-loop 
contributions in the two approaches in fact coincide. The remaining numerical difference 
between the diagrammatic calculation in the MS scheme and the result obtained by the 
EFT approach can thus be identified with new threshold effects due to non-logarithmic 
two-loop terms contained in the diagrammatic result. We furthermore show that in the 
analytic approximation employed in this paper, the dominant numerical contribution 
of these terms can be absorbed into an effective one-loop expression by choosing an 
appropriate scale for the running top-quark mass in different terms of the expression. 

This paper is organized as follows. To simplify the analysis, we focus completely 
on the radiatively corrected Higgs squared-mass in the "leading mf approximation" (in 
which only the dominant loop corrections proportional to mfhf ~ Gprn^ are kept). 
In addition, we choose a very simple form for the stop squared-mass matrix, which 
significantly simplifies the subsequent analysis while maintaining the most important 
features of the general result. In section 2 we sketch the derivation of the EFT result 
for the radiatively-corrected Higgs mass of refs. [p|-[TT[] at Olmfhfas) in the limit of 
mA ^ rriz- The corresponding result of the two- loop diagrammatic computation, under 
the same set of approximations, is outlined in section 3. In order to compare the two 
results, we must convert on-shell quantities to MS quantities. In section 4 we derive 
the relations between the on-shell and the MS values of the parameters m^, Xt and 
Ms in the limit of large Ms- Details of the exact calculation are given in Appendix A, 
while explicit relations between the on-shell and the MS parameters up to (9(m^/M|) 
are given in Appendix B. In section 5 the diagrammatic on-shell result is expressed 
in terms of MS parameters and compared to the result of the EFT computation of 
section 2. The logarithmic contributions are shown to coincide, and the remaining 
difference caused by non-logarithmic two-loop terms is analyzed. We argue that the 
remaining difference can be minimized by improving the EFT computation by taking 
into account the stop-mixing threshold contribution to the running top-quark mass. In 
addition, we demonstrate that in a simple analytic approximation to the Higgs mass, one 
can absorb the dominant two-loop contributions into an effective one-loop expression. 
In section 6, we summarize our results and discuss suggestions for future improvements. 
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2 Effective Field Theory Approach 



At the tree level, the mass matrix of the neutral CP-even Higgs bosons in the basis of 
weak eigenstates of definite hypercharge —1 and +1 respectively can be expressed in 
terms of mz, rriA and tan/5 = V2/vi as follows: 

^2,tree _ / sin^ /5 + m| cos^ /? -(m^ + m|) sin/?cos/? \ 

^ \^ — (m^ + m|) sin/9cos/? cos^ /? + m| sin^ /3 J' ^ ' 

Diagonalizing this mass matrix yields the tree-level prediction for the lightest neutral 
CP-even Higgs-boson mass 



2, tree 1 
= 2 



+ — J {m\ + m|)2 — Arrizm^ cos^ 2/3 



(2) 



For simplicity, we shall consider the limit of ^ mz and large supersymmetry- 
breaking masses characterized by a scale M5. Then, at energy scales below Ms, the 
effective low-energy theory consists of the Standard Model with one Higgs doublet. 
The corresponding Higgs squared-mass at tree-level is given by mf^^'^^^ = m\ cos^ 2/3. 
The dominant contributions to the radiatively-corrected Higgs mass within the EFT 
approach is based on the evaluation of the effective quartic Higgs self-coupling. A, eval- 
uated at the scale Q = frit- The value of X{Ms) is fixed by the supersymmetric boundary 
condition, although this value is slightly modified by one-loop threshold effects (denoted 
below by AthA), due to the decoupling of squarks at Ms with non-zero mixing. One 
then employs the Standard Model RGEs to obtain X{Q). Finally, the Higgs mass is 
obtained via m| = 2A(mt)u^(mt) [where v = 174 GeV is the Higgs vacuum expectation 
value] . The mass mj denotes the running top-quark mass in the MS scheme at the scale 
rrit. It is related to the on-shell (or pole) top-quark mass Mt = mf^ by the following 
relation ^ 

mt = mfi^imt) = 4 ' . , (3) 
i + -^cusyivit) 

where we have only included the QCD corrections to leading order in as- In eq. (^, 
the subscript 'SM' indicates that the running mass is defined in the usual way, i.e. in 
terms of the pure Standard Model (gluonic) contributions in the (modified) minimally 
subtracted dimensional regularization (DREG) scheme [ffH . 



The most important contributions to the mass of the lightest CP-even Higgs boson 
arise from the t-t sector of the MSSM, which is characterized by the following squared- 
mass matrix 

U2 _( Ml+mi + cos 2/5(1 - lsl,)ml m^X^ \ 

~ 1^ mtXt + ^? + I cos 2/5s2,m| ) ' 

where Xf = At ~ fi cot f3. The t-masses m^^ , m^^ the mixing angle 61 are determined 
at tree-level by diagonalizing A4^. Neglecting the numerically small contributions pro- 
portional to m| in the stop squared-mass matrix and setting 

M^^ = = MsusY, M| = M|usY + "^? (5) 
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leads to the simplified mass matrix 



In this approximation, the t-masses and the mixing angle are given by 

ml=Ml + \mM (7) 



i2 

f forXt<0 
f forXt>0, 



9: 



where by definition, m^^ < m^^. 

The one-loop threshold corrections to the quartic Higgs self-coupling, induced by 
the decoupling of stops, lead to a change of the effective quartic Higgs self-coupling at 
the scale Ms, 

X{Ms) = 1(9^ + cos^ 2/5 + AthA, (9) 

where the first term is the tree-level value of the quartic Higgs self-coupling in the 
effective low-energy Standard Model and the second term is the effect of the one-loop 
threshold corrections at the scale Ms |]7|,pil|]. 



A-A = i-,h\ \ [h\ - !(/ + (g) - ^,h\ (^) } + • • • , (10) 



where all couplings in eq. ( |10|) should be evaluated at the scale Ms- The running 
Higgs-top quark Yukawa coupling is related to the MS running top quark mass: 

mi(/i) = ht{y)v{[i) , (11) 

where the running Higgs vacuum expectation value, v^iMs) = v'^ijrit) i^'^ijrit), is gov- 
erned by the Higgs field anomalous dimension 



am) = 1 + 7;7^/^?(^t) In =1 • (12) 



^-^h^mt) In — 
In obtaining eq. (p!oD, an expansion in the variable 



1 1 

vfir — m 



= = ^ ^, < A,- < 1, (13) 

has been performed. Terms not explicitly exhibited in eq. (|TUp denote the contributions 
from higher powers in rrit/Ms and X^/Mg, which arise from the contributions of the t-t 
sector. Contributions from other supersymmetric-breaking sectors have been omitted 
for simplicity of the presentation. These contributions typically contribute no more 
than a few GeV to the radiatively-corrected Higgs mass. 
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As it was shown in refs. ||T0[ and [|Tl|], one can obtain the two-loop leading- logarithmic 
correction by expanding the parameter A up to order [ln(M|/m^)]^, 



A(mt) = XiMs) - (3xiMs)t + lP'^imt)t^ + ■ 



where \{Ms) is given by eq. @ and 
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m 



(14) 



(15) 



Following ref. |I0|, we define f3\ = ax\ + bx. Therefore 

X{mt) = X{Ms) [1 - ax{mt) t] - b^irfit) t [1 - ax{mt) t] - l(3'^{mt)t^ . (16) 
Here, 1 — axijnt) t = ^^"^(mf), where ^ is the Higgs field anomalous dimension [eq. (p!2D ]. 



Multiplying eq. (|T^) by 2t>^(mi), we obtain an equation for the Higgs squared- mass in 
the low-energy theory, which takes the following form: 



,2 /— 



(17) 



which defines the quantity A^d ^K^^i^t)- In eq. ([T^), 

mliMs) = 2X{Ms)v\Ms) , (18) 

where X{Ms) is given in eq. (^ with all couplings and masses evaluated at the scale Ms- 
In the present analysis, we are working in the approximation of hb = g = g' = 0. 
That is, we focus only on the Higgs-top quark Yukawa and QCD coupling effects. The 
relevant /?-f unctions for A, g^ and at scales below the scale Ms are given by [|19[ 



Wn^px = 6iX^ + Xhl - ht 
16n'p,2 ^ hi [Ihl - 8gl 
lQ7r'(3g2 = (-11 + iNf 



hi 

87r2 



15/if 



16^1 



9t. 



(19) 
(20) 
(21) 



where (3x = dX/ d In and Nf is the number of quark flavors with masses less than Q 
{e.g., Nf = 6 for scales between and Ms). Observe that we have included the domi- 
nant strong gauge coupling two-loop contribution to the j3 function of the quartic Higgs 
self-coupling, since it will contribute once we include all two-loop leading-logarithmic 
corrections. 

Using the above expressions, it is simple to flnd an approximate formula for the 
lightest CP-even Higgs mass in the large ttia limit. First, one obtains 



4 



47r^ v'^{mt) 



t 



1 + 



1 



1671' 



{Ihl - 32 



vra, 1 1 



^ — m? 



(22) 



To complete the 



where all couplings in eq. ( |22D are evaluated at the scale Q 
computation of m\{mt) [eq. (1T7|)], one must evaluate m\{Ms) [see eq. (|^)] in terms of 
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low-energy parameters. This is accomplished by using one-loop renormalization group 
evolution to relate X{Ms)v'^{Ms) to \{jnt)v'^{mt). In this way, one finally arrives at the 
expression]^ 



2, tree 



+ 



+ 



ml 



47r2 



t + 



167r2 



-^- 327r«, 



_xr 



2X1 
-Ml 



12M|^ 



+ 



4a, 



5/1? 



t . (23) 



The last two terms in eq. (|23|) refiect the two-loop single logarithmic dependence induced 
by the two-loop /3-function contribution to the running of the quartic Higgs self- coupling. 
It is interesting to note that these two terms are numerically close in size, and they tend 
to cancel each other in the computation of the Higgs mass. Eq. (PB]) differs from the one 
presented in ref. only in the inclusion of these terms, which although sub-dominant 
compared to the remaining terms, should be kept for comparison with the diagrammatic 
result. 

The full two-loop corrections to m| at 0{mfhj) have not yet been calculated in the 
diagrammatic approach; thus we neglect terms of this order in what follows.^ With a 
slight rewriting of eq. (|23|) we finally obtain the expression that will be compared with 
the diagrammatic result in the following sections: 



m 



2,tree , 3 Gpy/^ 4 



^ 2" TT- 



-3 



m. < In 




Ml 



(24) 



where we have introduced the notation Ms,Xt to emphasize that the corresponding 
quantities are MS parameters, which are evaluated at the scale /i = M5: 



Ms = Mf''iMs), Xt = X^%Ms 



MS/ 



(25) 



and irit 



m 



MS 
t,SM 



'rrit) as defined in eq. (||). 



3 Diagrammatic calculation 

In the diagrammatic approach the masses of the CP-even Higgs bosons are obtained by 
evaluating loop corrections to the h, H and /ii7-mixing propagators. The masses of the 
two CP-even Higgs bosons, nth and uihi are determined as the poles of this propagator 
matrix, which are given by the solution of 



2 2, tree 



^hH{q^ 



0, 



(26) 



^In the "leading approximation" that is employed here, there is no distinction between mhimt) 
and the on-shcU (or pole) Higgs mass, to/j. 

■^As noted below eq. (|2^), terms of O^m^h^) can be as numerically important as terms oiO{mfhfas). 
Hence, in a complete phenomenological analysis, one should not neglect terms of the former type. 
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where T^Hui.'f'), ^hH^'f) denote the renormahzed Higgs boson self-energies. In 

ref. |]13[ the dominant two-loop contributions to the masses of the CP-even Higgs bosons 
of 0{aas) have been evaluated. These corrections, obtained in the on-shell scheme, have 
been combined in refs. and ||T5[ with the complete one-loop on-shell result of ref. [ffl 



and the two-loop corrections of 0{mfh^) given in refs. p|-|TI|. 



The diagrammatic two-loop calculation of ref. [1^ involves a renormalization in the 
Higgs sector up to the two-loop level and a renormalization in the stop sector up to 
0{as). In the on-shell scheme, the renormalization in the stop sector is performed such 
that the t-masses m^^, m^^ correspond to the poles of the propagators, i.e. 

Re {ml ) = 0, Re S,-,,-, (m| ) = (27) 
for the renormahzed self-energies. In Ref. [|l^ the renormalization condition 



ReS,-^,-,K) = (2^ 



has been chosen to define the stop mixing angle.0 In ref. [|T6| a compact analytic 



approximation has been derived from the rather complicated diagrammatic two-loop 
result by performing an expansion in [eq. (|l^)] of the t-t sector contributions. 

The diagrammatic two-loop corrections to the Higgs mass also depend non-trivially 
on the gluino mass, which is a free input parameter of the supersymmetric model. In the 
EFT approach described in section 2, the gluino is decoupled at the same scale as the 
stops. Thus, in order to compare the results of the EFT and diagrammatic approaches, 
one must take rrig ~ 0{Ms). In this paper, we have chosen 



rrig = MsusY = V^l " '^t ■ (29) 

For the one-loop contributions from the other sectors of the MSSM the leading loga- 
rithmic approximation has been used [0, [111. In this approximation, the momentum 



dependence in eq. (P6|) is neglected everywhere. The resulting expression can thus be 
written as a correction to the tree-level mass matrix [eq. (|I|)]. The expression for ml in 
this approximation is obtained by diagonalizing the loop-corrected mass matrix. The 
compact analytic expression derived in this way, which is valid for arbitrary values of 
mA, has been shown to approximate the full diagrammatic result for mt rather well. 



typically within about 2 GeV for most parts of the MSSM parameter space . 

In the following we will restrict ourselves to the contribution of the t-t sector. In or- 
der to perform a simple comparison with the EFT approach of section 2, we only consider 
the dominant one-loop and two-loop terms of 0{m^h1) and 0{m'^h'^as), respectively. 
We focus on the case m^ ^ mz, for which the result for ml can be expressed in a 
particularly compact form. 



2 2,tree , 2,a , 2,a«s /on\ 



^In this paper, our analysis is presented in a simplified model of stop mixing, where the tree-level 
stop squared-mass matrix given by eq. (|^). In this case, ti and ^2 are states of definite parity at 
tree- level. Since parity is preserved to all orders in as, it follows that ^^^^^(p^) = (when electroweak 
corrections are neglected) and eq. (^ is trivially satisfied. 
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and neglect the non-leading terms of (9(m|/m^). Moreover, assuming that Ms S> Mt 
and neglecting the non-leading terms of 0{Mt/ Ms) and 0{m\/M^), one obtains the 
following simple result for the one-loop and two-loop contributions 



2,a 



nil 



71^ 



Mt 



In 



'Mf 
Ml 



+ 



XI 

Ml 




12 Ml 



-3 



a. 



MflW 



TT 



'Ml 
Ml 



2 + 



xl 

Ml 



In 



'Ml 
Ml 




(31) 



(32) 



The corresponding formulae, in which terms up to 0{Mf /Ml) are kept, can be found 
in Appendix B [see eqs. ( |B.1| ) and ( [B.2| )]. 

In eqs. (0) and (|32D the parameters Mt, Ms, Xt are on-shell quantities. Using 
eq. (^, the on-shell result for m| [eqs. (|30D-(|32D] can easily be rewritten in terms of the 
running top-quark mass Trif While this reparameterization does not change the form 
of the one-loop result, it induces an extra contribution at 0{aas)- Keeping again only 
terms that are not suppressed by powers of rrit/Ms, the resulting expressions read 



2,« 



3G^V2. 



-m. 



TC 

Xt 



In 



nit ' 



+ 



GFV^as 4 

-3 m. 



mi 



XI 

Ml 



+ 



4 

+ - - 2 

3 Ms 



8X1 17 Xf 
3M| ^ 36M| 




(33) 



(34) 



in accordance with the formulae given in ref. |T^ . 

We now compare the diagrammatic result expressed in terms of the parameters nit, 
Ms, Xt [eqs. (H) and (H)] with the EFT result [eq. (||)] which is given in terms of the 
MS parameters nit, Ms, Xt [eq. (^5])]. While the Xt-independent logarithmic terms 
are the same in both the diagrammatic and EFT results, the corresponding logarithmic 
terms at two-loops that are proportional to powers of Xt and Xt, respectively, are 
different. Furthermore, eq. (0) does not contain a logarithmic term proportional to 
Xt, while the corresponding term proportional to Xf appears in eq. (p^. To check 
whether these results are consistent, one must relate the on-shell and MS definitions of 
the parameters Ms and Xt. 

Finally, we note that the non-logarithmic terms contained in eq. ( ^4]) correspond to 
genuine two- loop contributions that are not present in the EFT result of eq. (^4|) . They 
can be interpreted as a two-loop finite threshold correction to the quartic Higgs self- 
coupling in the EFT approach. In particular, note that eq. (plj) contains a term that is 
linear in Xt This is the main source of the asymmetry in the two-loop corrected Higgs 
mass under Xt — *■ —Xt obtained by the diagrammatic method. The non- logarithmic 
terms in eq. (0) give rise to a numerically significant increase of the maximal value of 
nih of about 5 GeV in this approximation. 
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4 On-shell and MS definitions of Ms cind Xt 



Since the parameters p = {m?^, m?^, 6^, rrit} of the t-t sector are renormahzed differently 
in different schemes, the parameters Ms and Xt also have a different meaning in these 
schemes. In order to derive the relation between these parameters in the MS and in 
the on-shell scheme we start from the observation that at lowest order the parameters p 
are the same in both schemes, i.e. p = = p^^ in lowest order. Expressing the bare 
parameters in terms of the renormahzed parameters and the counterterms leads to 

pMs^^pMs^pOs^^pOS_ (35) 

Here 6p'~'^ is the on-shell counterterm in D = 4 — 2e dimensions, and according to 
the MS prescription Sp^^ is given just by the pole part of i.e. the contribution 

proportional to 1/e — 7e + ln47r, where 7e is Euler's constant. The MS parameters are 
thus related to the on-shell parameters by 

where Ap = 6p'~'^ — Sp^^ is finite in the limit D — 4 and contains the MS scale /i 
which can be chosen appropriately. In this paper, we only need to know Ap to 0{as) 
one-loop accuracy. In the following we will compare the result for expressed in terms 
of the on-shell parameters with results for rrih in terms of the corresponding MS 
parameters p^^, which are related to p*-*^ as in eq. (p^). 

In the EFT approach, the parameters Ms and Xt are running parameters evaluated 
at the scale fi = Ms [eq. (^)]. In the simplified model for the stop squared-mass matrix 
given by eq. @, the relations between the parameters Xt and Ms in the on-shell and 
MS scheme are obtained using 

r„2,OS ^ ^2,OS ^ ^^^OS ^ ^|OS ^ ^2,OS ^ ^^^OS ^ (37) 

mf^ = MIt mt{Ms)Xt , mf^ = M| ± mt{Ms)Xt , (38) 

where we have written Tnt{Ms) = mf^^{Ms) and Mt = m^^ as in section 2. In both 
eqs. (^) and (0), the upper and lower signs refer to > and X^^ < 0, respectively. 
In the model of stop mixing under consideration, there is no shift in the scalar top mixing 
angle to all orders in a^, from which it follows that |^p^| = \6^^\. 

Inserting the relation between m? , in the on-shell and the MS scheme into 
eqs. (^TD and (|38|) yields up to first order in as 

M| = M|'°^ + i (Am| + Am|) , (39) 

Xt = xO^^^Vt ± ^ (Am? - Ami ) , (40) 
* TrTt{Ms) 2mt V *2 hj ^ v ; 

where again the upper and lower sign in the last equation refers to X^^ > and 
X^^ < 0, respectively. In the second term of eq. (^0]) it is not necessary to distinguish 



(at one-loop) between mt{Ms) and M^, since Amj^, '^"^l ^i^is) quantities; hence, 
the generic symbol mt is used here. 
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In Appendix A, we have obtained explicit results for Am? , Am?^ and Mt/mt{Mi 



Inserting the appropriate expressions for these quantities into eq. (0), one observes that 
the functional form for Xt is the same for X^^ > and X^^ < [i.e., the sign difference 
in eq. (^) is compensated by the term (Am?^ — Am?^)]. As a result, it is no longer 
necessary to distinguish between these two cases. The case X^^ = (which formally 
would have to be treated separately) is understood as being included in eq. (|40|). 

Using the expansions given in Appendix A and setting the gluino mass according to 
eq. (p9|), we obtain to leading order in rrit/Ms 



Ml 



M. 



2, OS 



— -M^ 



X, = X, 



OS 



3 TT 

Mt 



S ) 



3 TT 



(41) 
(42) 



As previously noted, it is not necessary to specify the definition of the parameters that 
appear in the 0{as) terms. Thus, we use the generic symbol M| in the 0{as) terms 
of eqs. (p|)-(|4^). The corresponding results including terms up to O [mf/Mg) can be 
found in Appendix B. 

Finally, we need to evaluate the ratio Mt/fnt{Ms). The relevant expression is given 
in eq. ( |A.14| ). Using the expansions given at the end of Appendix A, we find to leading 
order in rrit/Ms 



mt{Ms) = nit 



1 + — In 



TT 



+ 



as Xt 
37r Ms 



(43) 



where rrit 



MS 



where terms up to O {mf / M| 



Tnpg-^[mt) is given in terms of Mt by eq. (|^). The corresponding formula, 
are kept, can be found at the end of Appendix B. Note 
that the term in eq. (|43| ) that is proportional to Xt is a threshold correction due to the 
supersymmetry-breaking stop- mixing effect. Inserting the result of eq. (^) into eq. 
yields: 

4Xt X2 



Xt 



X?' + ^Ms 



+ 



Ms 



Ml 



3Xt 
Ms 



In 



(44) 



It is interesting to note that Xt ^ when X^ = 0. Moreover, it is clear from eq. ( |4^ ) 
that the relation between Xt defined in the on-shell and the MS schemes includes a 
leading logarithmic effect, which has to be taken into account in a comparison of the 
leading logarithmic contributions in the EFT and the two-loop diagrammatic results. 

The above results are relevant for calculations in the full theory in which the effects 
of the supersymmetric particles are fully taken into account. However, in effective field 
theory below Ms, one must decouple the supersymmetric particles from the loops and 
compute with the Standard Model spectrum. Thus, it will be useful to define a running 
MS top-quark mass in the effective Standard Model, 
by: 



mfgjy[(/x), which to 0{as) is given 



m 



MS 
t,SM 



mt 



a.. 



1 + ^ln 



TT 



mr 



(45) 



At the scale Ms, we must match this result onto the expression for rnt{Ms) as computed 
in the full theory [eq. (1^)]. The matching is discontinuous at /i = Ms due to the 
threshold corrections arising from stop mixing effects. 
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In comparing with the EFT results of refs. |p|-[TT[l, it should be noted that the 
threshold correction in eq. (|43D were omitted. This is relevant, since TrTt{Ms) [or the 
related quantity ht{Ms), see eq. (p])] appears in the threshold correction to the quartic 
Higgs self-coupling X{Ms) [eqs. (^) and (|TD|)]. In refs. P-[TI|], mt{Ms) is re-expressed in 
terms of mt{mt) by using eq. (^Sf) rather than eq. (0). As a result, a two-loop non- 
logarithmic term proportional to Xt is missed in the computation of m/^. Such a term is 
of the same order as the two-loop threshold correction to the quartic Higgs self-coupling, 
which were also neglected in refs. |]8Hll|. However, in this work we do not neglect the 



latter. Hence, it would be incorrect to use eq. (|i5| ) in the evaluation of mt{Ms)- In 
section 5 we will apply rnt{^) with different choices of /i for the X^-independent and 
Xf-dependent contributions to m|, which will prove useful for absorbing numerically 
large two-loop contributions into an effective one-loop result. In the spirit of EFT, we 
will argue that for /i = M5, one should use the results of eq. (^Hf ) while for /x < M5, one 
should use eq. (|^). 

A remark on the regularization scheme is in order here. In effective field theory, 
the running top-quark mass at scales below Ms is the SM running coupling [eq. (^5])], 
which is calculated in dimensional regularization. This is matched onto the running 
top-quark mass as computed in the full supersymmetric theory. One could argue that 
the appropriate regularization scheme for the latter should be dimensional reduction 
(DRED) |^0[, which is usually applied in loop calculations in supersymmetry.]^ The 
result of such a change would be to modify slightly the two-loop non-logarithmic contri- 
bution to nih that is proportional to powers of Xt. Of course, the physical Higgs mass is 
independent of scheme. One is free to re-express eqs. ( pT]) and ([3^ ) [which depend on the 
on-shell parameters Mf, M5, Xt\ in terms of parameters defined in any other scheme. 
In this paper, we find MS-renormalization via DREG to be the most convenient scheme 
for the comparison of the diagrammatic and EFT results for mh. 



5 Comparing the EFT and diagrammatic results 



In order to directly compare the two-loop diagrammatic and EFT results, we must 
convert from on-shell to MS parameters. Inserting eqs. (|4l|) and (|4^) into eqs. (^3D and 
fBl), one finds 



2,0 



3 GfV2_ 



2, acts 




(46) 



(47) 



Comparing eq. (47) with eq. (0) shows that the logarithmic contributions of the 
diagrammatic result expressed in terms of the MS parameters m^, Ms, Xt agree with 



^In order to obtain the corresponding DRED result, one simply has to replace the term Aas/Sn in 
the denominator of eq. (0) by bas/Sn. 
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the logarithmic contributions obtained by the EFT approach. The differences in the log- 
arithmic terms observed in the comparison of eqs. (^) and (|3^ ) with eq. (^) have thus 
been traced to the different renormalization schemes applied in the respective calcula- 
tions. The fact that the logarithmic contributions obtained within the two approaches 
agree after a proper rewriting of the parameters of the stop sector is an important con- 
sistency check of the calculations. In addition to the logarithmic contributions, eq. (47) 
also contains non-logarithmic contributions, which are numerically sizable. 

In fig. [^, we compare the diagrammatic result for ruh in the leading approximation 
to the results obtained in section 2 by EFT techniques, for two different values of tan j3. 
However, as noted at the end of section 4, in the derivation of the EFT result of eq. (|2^ ) 
the supersymmetric threshold corrections to mt{Ms) were neglected. Thus, Xt, which 
appears in eqs. (|^) and (^T]) , is not precisely the same as the Xt parameter appearing 
in the EFT result of eq. due to the difference in the definition of mt{Ms) [eq. 
and mff-^ [eq. (^Sf)]. Taking this difference into account in eq. (|42|) , it follows that the 
Xf parameter that appears in eq. ( p4D is given by X/ = [1 + {as/3TT){Xt/Ms)]. It 
can be easily checked that the change from Xf to X^ does not affect the comparison 
of two-loop logarithmic terms between eqs. (p^) and (^Tj). Moreover, the difference 
between Xt and Xj. is numerically small. In fig. |I|, the diagrammatic result for rrih is 
plotted versus Xt, while the EFT result is plotted versus Xj.. 

While the diagrammatic result expressed in terms of mt, Ms, Xt agrees well with 
the EFT result in the region of no mixing in the stop sector, sizable deviations occur for 
large mixing. In particular, the non-logarithmic contributions give rise to an asymmetry 
under the change of sign of the parameter Xt, while the EFT result is symmetric under 
Xt — >■ —Xt- In the approximation considered here, the maximal value for in the 
diagrammatic result lies about 3 GeV higher than the maximal value of the EFT result 
for tan/3 = 1.6. The differences are slightly smaller for tan/? = 30. In addition, as 
previously noted, the maximal-mixing point {Xt)ma,x [where the radiatively corrected 
value of rrih is maximal] is equal to its one-loop value, {Xt)ma.x — ^VQMs, in the EFT 
result of eq. (p^ , while it is shifted in the two-loop diagrammatic result. However, fig. |l] 
illustrates that the shift in {Xt)ma.x from its one- loop value, while significant in the two- 
loop on-shell diagrammatic result, is largely diminished when the latter is re-expressed 
in terms of MS parameters. 

The differences between the diagrammatic and EFT results shown in fig. || can 
be attributed to non-negligible non-logarithmic terms proportional to powers of Xt. 
Clearly, the EFT technique can be improved to incorporate these terms. As previously 
discussed, one can account for such terms in the EFT approach by: (i) including one- 
loop finite threshold effects in the definition of (M5), and (ii) including two-loop finite 
threshold effects to AthA [eq. (jl^)] in the matching condition for A(M5).| Although both 
effects are nominally of the same order, it is interesting to investigate what fraction of 
the terms proportional to powers of Xt in eq. (|47|) can be obtained by including the 
threshold effects in the definition of rnt{Ms). 

^An additional two-loop 0(/ijas) correction to X{Ms) can arise because the Higgs self-coupling in 
the MS scheme does not precisely satisfy the supersymmetric relation [X{Ms) — \{g^ + 5'^) cos^ 2/3] in 
the supersymmetric limit. This correction corresponds to a matching of the MS and DR couplings at 
the scale M5 We will count this as part of the two-loop finite threshold effects. 
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Leading m,'' approximation at 0(a a^) 
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Figure 1. Comparison of the diagrammatic two-loop 0{m^hfas) result for rrih, to leading order in 
mt/Us [eqs. (H) and @] with the EFT result of eq. (||). Note that the latter omits the one-loop 
threshold corrections due to stop mixing in the evaluation of rnt{Ms)- Since this quantity enters in the 
definition of Xt [see eq. (|4^)], the meaning of Xf^^ plotted along the a:-axis is slightly different for the 



diagrammatic curve, where Xf^^ = Xt, and the EFT curve, where X}^^ ^ Xt[l + (as/37r)(Xt/Ms)]. 
See text for further details. The two graphs above are plotted for Ms — tua — (m| -)-TO()^/^ — 1 TeV 
for the cases of tan (3 = 1.6 and tan/3 = 30, respectively. 
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In ref. it was shown that the leading two- loop contributions to m\ given by the 
EFT result of eq. (p^ ) could be absorbed into an effective one-loop expression. This was 
accomplished by considering separately the X^-independent leading double logarithmic 
term (the "no-mixing" contribution) and the leading single logarithmic term that is 
proportional to powers of Xt (the "mixing" contribution) at 0{mfhfas). Both terms 
can be reproduced by an effective one-loop expression, where mt in eq. (^6|), which 
appears in the no-mixing and mixing contributions, is replaced by the running top- 
quark mass evaluated at the scales fit and /i^, respectively: 

no mixing: jjt = {7ntM s)^^"^ , mixing: ^i = Ms- (48) 

That is, at 0{mfhf as), the leading double logarithmic term is precisely reproduced 
by the single-logarithmic term at Olmfhf), by replacing with mt{fit), while the 
leading single logarithmic term at two-loops proportional to powers of Xt is precisely 
reproduced by the corresponding non-logarithmic terms proportional to Xt at 0{mfhf), 
by replacing rrit with mt{Ms)- 

Applying the same procedure to eq. ( ^7|) and rewriting it in terms of the running 



top-quark mass at the corresponding scales as specified in eq. (|i8|), we obtain 



= In (^) + mf(M.)§ ( 1 - ^ U , (49) 



■5 / ^"S 

7r2 TT 1 6 I Mg ] + I 3 9 Mg ' 36 



= -3^^m^ ^ ( ^ ) + ^ ( ^ - + ) ) . (50) 




Indeed, the Xj-independent leading double logarithmic term and the leading single log- 
arithmic term that is proportional to powers of Xt have disappeared from the two-loop 
expression [eq. (|5D|)], having been absorbed into an effective one-loop result [eq. (^5])] 
(denoted henceforth as the "mixed-scale" one- loop EFT result). Of the terms that re- 
main [eq. (^np], there is a subleading one-loop logarithm at two-loops which is a remnant 
of the no-mixing contribution. But, note that the magnitude of the coefficient (1/6) 
has been reduced from the corresponding coefficients that appear in eqs. (^2|) and (|47| ) 
[—2 and 2/3, respectively]. In addition, the remaining leftover two- loop non- logarithmic 
terms are also numerically insignificant. We conclude that the "mixed- scale" one- loop 
EFT result provides a very good approximation to m\, in which the most significant 
two-loop terms have been absorbed into an effective one-loop expression. 

To illustrate this result, we compare in fig. ^ the diagrammatic two- loop result 



expressed in terms of MS parameters [eqs. (^) and (^Tf)] with the "mixed-scale" one- 



loop EFT result [eq. (^)] as a function of Xt. In order to make a fair comparison of 



two-loop expressions, we first evaluate eq. (^91) as a perturbation expansion which is 
truncated beyond the 0{a1) term.Q It is this result that is plotted as a dashed line in 
fig. ||. Note that by construction, the sum of the two-loop truncated version of eq. (|49| ) 
and the leftover two-loop term given by eq. (^) is equal to the sum of eqs. (^61) and 
(^). That is, the difference between the solid and dashed lines of fig. ^ is precisely 

"^For example, using eq. (^), we would write (/ij) = [l + 4(as/7r) ln(mt /yu^ )] and insert this 
into eq. 
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Leading m,'' approximation at 0(a a^) 





Figure 2. Comparison of the diagrammatic two-loop 0{mjhfas) result for rrih, to leading order 
in rnt/Ms [eqs. (H) and @] with the "mixed-scale" one-loop EFT result [eq. (|9|)]. Note that 
the latter now includes the threshold corrections due to stop mixing in the evaluation of rnt{Ms) in 
contrast to the EFT results depicted in fig. |^. "Mixed-scale" indicates that in the no-mixing and 
mixing contributions to the one- loop Higgs mass, the running top quark mass is evaluated at different 
scales according to eq. (^). See text for further details. The two graphs above are plotted for 
Ms = rriA = {m- + rn^y/^ = 1 TeV for the cases of tan (3—1.6 and tan/3 — 30, respectively. 
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equal to the leftover two-loop term given by eq. (pOD , which is seen to be numerically 
small. Hence, within the simplifying framework under consideration [i.e., only leading 
t-t sector-contributions are taken into account assuming a simplified stop squared-mass 
matrix [eq. d)], with Ms, ruA ^ frit and rrig = Msusy), we see that the "mixed-scale" 
one-loop result for ruh provides a very good approximation to a more complete two-loop 
result for all values of In the EFT picture this means that, once re-expressed in 
terms of the appropriate MS running parameters, the dominant contributions to the 
lightest CP-even Higgs mass arising from the two-loop threshold corrections induced by 
the decoupling of the stops, have their origin in the one-loop threshold corrections to 
the Higgs-top quark Yukawa coupling. 

In the present analysis we have focused on the leading contributions of the t-t sec- 
tor of the MSSM. However, these contributions alone are not sufficient to provide an 
accurate determination of the Higgs mass (and can be off by 5 GeV or more in certain 
regions of the MSSM parameter space). In any realistic phenomenological analysis of 
the properties of the Higgs sector, one must include sub-leading contributions of the 
t-t sector as well as contributions from other particle/superpartner sectors. Such con- 
tributions have been obtained within the EFT approach in refs. [|10| and which 
incorporates one-loop leading logarithmic terms from all partner/superpartner sectors, 
plus single and double logarithmic two-loop contributions from the t-t and b-b sectors. 
The full one-loop diagrammatic result is known , and this along with the diagram- 
matic two-loop result from the t-t sector at 0{aas) are included in the Fortran code 
FeynHiggs |]22|. The impact of the additional contributions to the radiatively-corrected 
Higgs mass on phenomenological studies have been investigated in refs. |15,p3 . 



In the large tan (3 regime, the b-b sector is especially important. Here the corrections 
induced by the bottom Yukawa coupling become relevant, and one should correspond- 
ingly include the bottom mass corrections originating from the decoupling of the super- 
symmetric particles. These corrections are enhanced by a large tan (3 factor and hence 
can have a sizable impact on the phenomenology of the Higgs sector. Some of the most 



relevant consequences of these corrections have been recently discussed in refs. p4|-[26 



6 Discussion and Conclusions 

In this work, we have compared the results for the lightest CP-even Higgs-boson mass 
obtained from the two-loop 0{aas) diagrammatic calculation in the on-shell scheme 
with the results of an effective field theory approach. In the latter, the two-loop 0{aas) 
terms are generated via renormalization group running of the effective low-energy pa- 
rameters from the supersjTiimetry-breaking scale, Ms to the scale mj. We have focused 
on the leading 0{mfhfas) two-loop contributions to m| in the limit of large ttla and 
Ms- In this case, the effective field theory below Ms is the one-Higgs-doublet Standard 
Model, which greatly simplifies the calculation. In addition, the gluino mass was set to 
TUg = Msusy = {Mg — m^Y^'^. The resulting transparent analytic expressions for the 

^Strictly speaking, the analytic approximations of this paper break down when WitXt ~ M^. Thus, 
one does not expect an accurate result for the corresponding formulae when Xt is too large ||9|,[Tl|,[l6|. 
In practice, one should not trust the accuracy of the analytic formulae once Xt > {Xt)ina.^. 
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radiatively-corrected Higgs mass were well suited for investigating the basic relations 
between the various approaches. In order to compare the on-shell diagrammatic and 
effective field theory approaches, one must note two important facts. First, the two 
calculations are performed in different renormalization schemes. Hence, the resulting 
expressions actually depend on soft-supersymmetry-breaking parameters whose defini- 
tions differ at the one-loop level. Second, the diagrammatic calculation includes genuine 
non-logarithmic two-loop corrections to the lightest CP-even Higgs-boson mass. In the 
effective field theory approach, these would correspond to two-loop threshold corrections 
resulting from the decoupling of the two heavy top squarks in the low energy effective 
theory. 

Previous comparisons of the corresponding two-loop results revealed an apparent dis- 
crepancy between terms that depend logarithmically on M5, and different dependences 
of the non-logarithmic terms on Xt. However, after re-expressing the one- loop and the 
two-loop terms of the on-shell diagrammatic result in terms of MS parameters, i.e., ap- 
plying the same renormalization scheme for both approaches, we have shown that the 
discrepancy in the logarithmic dependence on Ms of both expressions disappears. This 
constitutes an important consistency check of the calculations. There remain, however, 
genuine non-logarithmic two-loop contributions in the diagrammatic result. They give 
rise to an asymmetry under Xt —Xt., while the effective field theory computations 
that neglected the two-loop threshold corrections due to stop mixing only yield results 
that are symmetric under the change of sign of Xt. Moreover, the non- logarithmic 
two-loop contributions of the on-shell diagrammatic computation, in the approxima- 
tions considered in this paper, can increase the predicted value of mh by as much as 
3 GeV.0 Finally, they induce a shift of the value of Xt where mh is maximal relative 
to the corresponding one-loop value (Xt)max — ±-\/6M5. It is interesting to note that 
this shift is more (less) pronounced when mu is expressed in terms of Xt in the on-shell 
(MS) scheme. The effect of the leading non- logarithmic two-loop contributions can be 
taken into account in the effective field theory method by incorporating the OQilag) 
Xt-dependent corrections into the boundary conditions of the effective quartic Higgs 
self-coupling at the scale Ms and performing a proper one- loop 0{as) matching of the 
running Higgs-top Yukawa coupling at M5. 

In ref. [|ll|, it was shown that the leading two-loop contributions to m\ could be 
absorbed into an effective one-loop expression by the following procedure. The run- 
ning top-quark mass that appears in the X^-independent one-loop expression for is 
evaluated at the scale Ht = {MsmtY/"^ . In the corresponding X^-dependent terms, the 
running top quark mass is evaluated at the scale ni = Ms- The result, which we call the 
mixed-scale one-loop EFT expression neatly incorporates the leading two-loop effects. 
In this paper, we have extended this result by explicitly including stop mixing effects in 
evaluating the running parameters fnt{Ms) and ht{Ms). By doing so, we are able to in- 
corporate some portion of the genuine leading non-logarithmic two-loop contributions. 
The remaining terms at this order would then be identified with two-loop threshold 
corrections to the effective Higgs quartic coupling at the scale Ms. Remarkably, the 
latter turn out to be numerically small. This means that the mixed-scale one-loop EFT 



^If only the on-shell top-quark mass is re-expressed in terms of the corresponding MS parameter, 
the resulting increase in rrih can be as much as 5 GeV. 
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expression for rrih provides a rather accurate estimate of the radiatively-corrected mass 
of the hghtest CP-even Higgs boson of the MSSM. 

The above results have been obtained in a rather simple setting. A special choice 
for the stop squared-mass matrix was made [eq. (|])] to simplify our analysis. The 
gluino mass was fixed to a value of order M5. The leading 0{mfhf) corrections were 
neglected. Subleading terms of 0(m|/ijas) and 0{m\h^) terms were also neglected. 
For example, consider the effect of varying the gluino mass. The two-loop diagrammatic 
results of refs. |lT3Hl6| showed that the value of irih changed by as much as ±2 GeV as 



a function of rrig, for rrit < rrig < Ms- The gluino mass dependence can be treated 
in the EFT approach as follows. Let us assume that Ms characterizes the scale of the 
squark masses, and rrig < Ms- Then, at scales below Ms one integrates out the squarks 
but keeps the gluino as part of the low-energy effective theory. However, since the 
gluino always appears with squarks in diagrams contributing to m| at two-loops, once 
the squarks are integrated out, they no longer affect the running of any of the relevant 
low-energy parameters below Ms- However, the gluino mass does affect the value of 
ht{Ms) and mt{Ms) [the relevant formulae are given in Appendix A]. Thus, in the EFT 
approach, gluino mass dependence enters via the threshold corrections to the Higgs-top 
quark Yukawa couplings. 

The case of a more general stop squared-mass matrix can be treated using the 
same methods outlined in Appendix A.[^ Here the computations are more complicated 
since there is now one- loop mixing between ti and ^2- In the EFT approach, one must 
decouple separately the two stops, and include the most general stop-mixing effects in 
the determination of the relation between the on-shell and MS parameters. The leading 
Oim^hf) corrections in the EFT approach can be incorporated as in section 2 [Ty,|TT]] by 



extending the computations of Appendix A to include the one-loop 0{h^) corrections 
to the running top-quark mass and stop sector parameters. However, at present, one 
cannot check these results against an 0{mfh^) two-loop diagrammatic computation, 
since the latter does not yet appear in the literature in full generality. 

Going beyond the approximations made in this paper, the next step is to incorporate 
the above improvements, as well as the next subleading contributions of 0{m%h1as) and 
0{m\hf) into the computation of m\. One might then hope to show that a complete 
mixed-scale one- loop EFT result, suitably generalized, provides a very good approxi- 
mation to the radiatively-corrected CP-even Higgs mass of the MSSM. Such an analysis 
could be used to organize the most significant non-leading one-loop and two-loop con- 
tributions to m\ and provide some insight regarding the magnitude of the unknown 
higher-order corrections, thus reducing the theoretical uncertainty in the prediction for 
rrth. This would have a significant impact on the physics of the lightest CP-even Higgs 
boson at LEP2, the upgraded Tevatron and the LHC 



^"The transformation from on-shcU to MS input parameters for the case of the most general stop 
squared-mass matrix will be included in the new version of the program FeynHiggs. 
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Note Added 



After this work was completed, we received a paper [^] which discusses many of the 
same issues that we address in this work. In ref. the two-loop effective potential 
of the MSSM is employed, including renormalization group resummation of logarith- 
mic terms, and the leading non-logarithmic two-loop terms of 0{aas)- The relation 
between on-shell and MS parameters are also taken into account. The end results are 
qualitatively similar to the ones obtained here. 
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Appendix A: Relations between on-shell and MS def- 
initions of mt, Ms cind Xt 

In this appendix, we derive the relations between the on-shell and MS definitions of 
rrit, Ms and X^. We have checked that these results agree with similar results given 
in ref. [Q. These results are derived in a model where the stop mass-squared matrix is 
given by eq. (|^) . The corresponding stop squared-masses and mixing angle are given by 
eqs. (|^ and (||). Note that in this model, the top-squark mass eigenstates, ti and t2, are 
states of definite parity in their interactions with gluons and gluinos. The corresponding 
Feynman rules are shown in fig. ^. 

Consider first the one-loop contribution to the top-quark two-point function at 0{as) 
due to: (i) the top-quark/gluon loop [fig. §(a)] and (ii) the stop/gluino loop [fig. ||(b)]. 
Divergences are regulated by dimensional regularization in D = 4 — 2e dimensions and 
removed by minimal subtraction. Including the tree-level contribution (which is equal 
to the negative of the inverse tree- level propagator) , the end result isQ 

r(2) (p) = ^y - - !^ {^[1 + 2B,{p'- ml 0)] - 2m,[l - 2So(/; m^, 0)] 

+ ^[Bi{p'; ml Ml - m^X,) + B^{p^; m?. Ml + m^Xt)] 
-mg[8,{plmlMl-mtXt)-Bi{plmlMl + mtXt)]} , (A.l) 

where T°'T°' = Cpl is the SU(3) quadratic Casimir operator in the fundamental repre- 

-'^^All results in this section are given in the MS subtraction scheme. In the DR scheme, aU the 
formulae of this appendix still apply except in the case of the top-quark gluon loop. To obtain the 



corresponding DR result, simply remove the additive factors of 1 in the two occurences in eq. (A.l). 
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Figure 3. Fcynman rules for top-squark interactions in a model where the stop mass-squared matrix 
is given by eq. (^. 
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Figure 4. One-loop contributions to the top quark mass. 



sentation, Cp = 4/3, and 



-1 



dyy"- In 



'mly + (1 - y) - - yY 



(A.2) 



The i3„ (n 



one-loop computations |28 



0,1) are related to the standard two-point loop functions that arise in 



2J 5 



2) ■> 



(A.3) 



where all occurrences of A = (47r)'^r(e) are removed in the minimal subtraction pro- 
cedure. Note that BQ{p'^]m\,m'^ is invariant under the interchange of m\ and m|, 
whereas 



Bi{p ]m^,mY) = -Bi{p ]m-^,m^) - Bq{p ]m-^,m^) . 



(A.4) 



In eq. ( [A.l| ), ^ is the arbitrary mass parameter of the MS-scheme. The on-shell top- 
quark mass, Mt, is defined by T^^\^ = Mt) = 0. It follows that 



Mt = mt{n) + 



An 



+ Bi{ml, ml M| - nitXt) + Bi{ml, m^, M| + m^X^ 



rrin 



Bo{ml, ml Ml - mtXt) - Boim]- ml Ml + m^X, 



(A.5) 



In the 0{as) terms above, we simply use the generic notation mt for the top-quark 
mass, since to one- loop accuracy one need not distinguish between agMt and asmt{^). 
As previously noted, the relation between the top-quark mass in the on-shell and DR 
schemes is obtained by dropping the —1 (which does not multiply a loop-function) in 
eq. (|A~5|) . 

Two of the loop functions are easily evaluated: BQ{mf;ml 0) = 2 - In(m2//i2) 
and Bi{mf;mlO) = — 1[3 — ln(m^//i^)]. We note the following curious fact. In the 
supersymmetric limit (m^^ = m^^ = and = 0) at one-loop, the relation between 
the on-shell and DR running top-quark mass evaluated at fi = mt is precisely the same 
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Figure 5. One-loop contributions to the top squark mass. 
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as the corresponding relation between the on-shell and MS running top-quark mass in 
non-supersymmetric QCD. 

Next, we examine the one-loop contributions to the top-squark two-point function 
at 0{as)- The contributing graphs are shown in fig. ^. We immediately note that graph 
(c) of fig. ^ vanishes in dimensional regularization. Moreover, since ti and ^2 are states 
of definite parity in a model where the stop mass-squared matrix is given by eq. (P), 
the one-loop mixing of ti and ^2 vanishes to all orders in a^- Including the tree-level 
contribution, the final result for the top squark two-point function isf^ 



nir 



IT 



for j = 1, 2, where 



, / m 



(A.6) 
(A.7) 



is related to the standard one-point loop function Ao{m'^) = m^A -|- Aolrn^). 

There is no distinction in this calculation between the MS and DR schemes. The 
on-shell stop squared-masses are defined by i'fjip'^ = {fn'?^)'^) = 0. Noting the form for 
the tree-level squared-masses [eq. (0)], it follows that: 



CpOis 



f{MlTmtXt)±g{MlTmXt 



where Mf is the on-shell top quark mass and 



f{p ) = Aoim-g) +mtBo{p ■,mt,m~g)+p Bi{p ;mi,m^ 
g{p'^) = m-gmtBoip'^; , m?) . 



(A.9) 



2 OS 

It is then straightforward to solve for Mg and MfX^ in terms of the corresponding 
MS quantities evaluated at the scale fi = Ms- Using the notation of eq. (|25D, we obtain 



2, OS 



M. 



CpCXs 



2-K 



f{Ml + mtXt) + f{Ml-mtX, 



-g{Ml + mtXt) + g{Ml - rritXt 



mt{Ms)Xt + 



CpOs 
r2 



fiMl + mtXt)-fiMl-mtXt) 



-g{Ml + mtXt) -g{M's- mtXt 



(A.IO) 



(A.ii: 



where all loop functions in eqs. ( |A.1(]| ) and ( |A.11| ) are evaluated at /i = Ms- Dividing 
eq. ( |A.11| ) by Mj, and using eq. ( [A.5| ) to evaluate mt{Ms) / Mt, one obtains a direct 
relation between X^^ and Xt- 

^^In deriving eq. ( A.6 ), we note that the contribution of fig. ||(d) [only the case i = j yields a non-zero 
contribution, which is equal to (zCi?as/47r)^o(™t )] cancels a similar term that arises in the evaluation 
of fig. 1(b). 
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(A.12) 



To obtain the expansions derived in Appendix B, we consider the case of m? 
MgusY = ~ ^t- introduce the following notation: 

Ms' '-Ms- 

We are interested in the limit of z ^ 1 and Xt < 1. First, consider the relation 
between the on-shell and running top quark mass. Using eq. (|A.2|) , we must evaluate 
the following integrals: 

rfyy" In [l ± xtzy - z\l - , (A.13) 



for n = 0, 1. Eq. ( [A.5| ) then yields 



mt{Ms 



Mt{l + 



Ait 



4 + 61nz - Ji 



{+) 



Jl ^ + -(1 — ^ ) yoQ 

(A.14) 

Expanding out the logarithm in the integrand of J^^^ in a double power series in xt and 
z and integrating term by term, one readily obtains the result given in eq. ( p.5| ). 

Second, consider the relation between the on-shell and MS definitions of Ms and 
Xt obtained in eq. ( |A.1CI|) . Using the integral expressions for the loop functions that 
appear in eq. ( |A.9| ), one must evaluate the following integrals: 



1 



dyy"" In y^{l - zxt) + yz{xt - 2z) + z^ 



(A.15) 



for = 0, 1. In this case, one cannot simply expand the logarithms about = z = 0, 
since the integration range extends down to y = 0. Instead, we have used Mathematica 
to evaluate the integral exactly, and then perform the double expansion in Xt and z. 
The result for Iq up to Oi^x'jz'^) is 



-2 + TTZ + z^{2\nz - 1) - \tiz^ + \z^ 



+Xt 



i~4 

z\nz + TTz^ + z^{2hiz - \] 



^-KZ 



+xl \-\t^z - z^{\\iz + 1) + ^Tiz^ + 2^(2 In z 

+xl [j^Z- \t^Z^ - Z^{hlZ+^) + ifTTZ 

-J-T,z + ^^z^-§.TTz'-z\\nz + %] 



(A.16) 



128 ' 12 256 ^ ' 12' 

We have checked the validity of this expansion using numerical integration. One can 
derive a similar expression for Ii either directly, or by noting that: 

1 



2(1 - zxt] 



z\l - 2 In ^) - (1 - z')[l - ln(l - z^)] - z{xt - 2z)Iq 



\- z^ihiz+l) + TTZ^ + z^{2 In z 



2 

+ Xt 

+ xl 

+ xl 
+ xt 



^TTZ 



z%3\nz + 2) + \ttz^ 



\z^{^^ z + 2)- fTT^^ - z\^ In ;z + f ) 

^-KZ^ + Z^{\nZ+^) - f TTZ 

-^^^ + ^7r/ + /(|lnz + f) 



(A. 17) 
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Inserting these results into eqs. ( A. 10 ) and ( A.ll|) , and expanding out the remaining 
factors [e.g., mgirit = Mgz{l — z'^Y^'^, etc.], one ends up with the results given in 
eqs. (il) and (plj). 



Appendix B: Results up to O {rn^/Mg 

We list the necessary formulae in order to derive the results given in sections 3 and 4 
up to terms of O (m^/Mg). As before, we define xt = Xt/Ms and z = rrit/Ms. In the 
approximation discussed in section we obtain expansions that are valid in the limit 
of z -C 1 and Xt < 1. 

The diagrammatic result in the on-shell scheme for the one-loop and two-loop con- 
tributions to ml, in the approximations used in this paper, is given up to 0{z*) by 

3GfV2 

1// ' — / in ^ -t- I — — ^~ ' 2 



m 



2,a 



2, acts 



2 7r2 



1 - ^x^ - z^(^ - ^x^) 



14 2 



} 



(B.i; 



M^iAln^z + 



-4 - 2x'f + \z^ (6 + ^Ixt + 'i'ixl - 26 

Inz 



X, 



+ (2 - 41xt - 10x^2 + 91Xi^ + 454 

- \xt{% - x\) + ^-nzxt (48 + 24xt - \^x\ - 7x 
+ ^60 — 6x4 — 106x^ — 6I2 



Y^nz^Xt 



(l92 + 72xt - 236x? 



^z^xt (57 + 34xt + 95xt + 20x^ 



lllxf) 
(B.2) 



Eq. ( [B.2| ) is a generalization of the corresponding formula given in ref. [^]. 

The relations between the MS parameters M5 and Xt [eq. (pS])] and the correspond- 
ing on-shell parameters M^^, X^^ can be obtained from eqs. ( |A.1U| ) and ( |A.11| ) using 
the expansions of eqs. ( |A.16| ) and ( [A.17| ). The end result up to 0{z^) is given by 



M| = M; 



2,0s 



2a. 
37r 



-4 - Y^^^ (12 + 24:Xt + 6x 



9™d _ ™4 



+ i7rz2xt(16 + Sxt - 2x^ - xf) + ^z^{6 - 6x? - 23x^ - lOx^ 



+ 



z^{2 - 2xt - xf) + z^xt{5 + 2xt - 2x1 - ^ 



In z 



(B.3) 



OS 



m 



OS 



2a 



+ -^Ms{ 4 - ^7r2(128 + Uxt - 16x? - Sxf - xf) 



37r 



+ 1^^X4(6 + 12xt + 4x^ - x^^) + ^7rz^(128 + 32xt - 128x^2 - 60x^ + 17xf) 



+ 



^z^(30 + Qxt - 6x^ - 23xf) 
2^(2 + Xt){-2 + xl) + z^{2 - hxl - 2x1 + 2a; 



In z 



(B.4) 



To complete the evaluation of Xt we need to find a relation between the on-shell 
top-quark mass, Mt, and the running top-quark mass, mt{Ms), evaluated at the scale 
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Ms- Using eq. ( A. 14 ), one obtains the following expansion: 



( (y. V / 
mt{Ms) = Mt\l + ^[-A + Q\nz + Xt + z'[\ + 1x1 + 1 

,74/l_J_ I 1^2 I J_ 3 I J_ 4 I J_ 5\1 \ 
~'~ VS 24 * ~'~ 6 < ~^ 12 t "T 15 * )\ / ' 



(B.5) 



Eq. ( p.5|) was derived using DREG. In order to obtain the corresponding formula using 
DRED (which yields a formula for the top-quark mass in the DR scheme), simply replace 
—4 with —5 in the first term of eq. (|B.5|) after the left square bracket. 

Note that eq. ( p.5|) provides a connection between rnt{Ms) and the on-shell mass 
Mt in the full supersymmetric theory. In the limit of large Ms with fixed Xt/Ms-, 
the threshold correction arising from the stop mixing effects does not vanish. On the 
other hand, the MS top-quark mass. Wit = ^Timi^t) is defined in the low-energy (non- 



supersymmetric) effective theory via eq. (y). Thus, in eq. (|B.5|) , we may replace Mt 
with rfit simply by removing the factor of — 4a;s/37r. To leading order in mt/Ms-, one 
immediately obtains eq. (^31). 
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